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The thermodynamic approach to non-equilibrium dynam- 
ics describes the state of macroscopic systems by means of 
a collection of intensities or intensive variables. The latter 
are by definition the differentials of the entropy with respect 
to the set of extensive constraints. The environment is di- 
rectly involved in controlling the intensities. The isolation 
paradigm is negated. The general principles substanciating 
the approach are restated and expanded to multi-component 
systems. The procedure is applied to the prediction of trans- 
port processes (viscosity and thermal conductivity) of mix- 
tures of atomic gases. Theoretical results are compared with 
published experimental data. 



I. INTRODUCTION 

Excepting recent work centred on trajectory calcula- 
tions jiJD, the common parenthood of most theoretical 
approaches to the evolution of irreversible phenomena 
traces back to Boltzmann's dynamic equation |^ . The 
core of the model is based on the isolation paradigm. In- 
teraction with the surroundings is claimed not to be the 
necessary condition for breaking the correlations repre- 
senting particular non-equilibrium initial conditions. By 
contrast, the molecular chaos hypothesis is presumed. 
Although this has been violently criticized in the early 
days, there is presently a considerable literature in math- 
ematical physics arguing in favour of this initial assump- 
tion Some claim that chaotic Hamiltonian dy- 
namics is the key to decorrelation leading to the required 
molecular chaos conditions 0. 

Contrasting with the former, some other schools of 
thoughts suggest that correlation breaking is a property 
typically imported from the surroundings |§HlCf|. They 
confirm that strictly isolated systems would follow Liou- 
villian dynamics with conservation of entropy. The ap- 
parent success of Boltzmann's equation and its different 
approximations is claimed to result from its a priori in- 
troduction of the molecular chaos assumption. The latter 
masks or hides non-Hamiltonian dissipativc components 
equally present in the global dynamics, as the effect of 
unavoidable interactions with the external world p!l|^2| . 

One of the arguments of the followers of the isolation 
paradigm is that transport properties characterizing fluid 
systems and typical of irreversible dynamics (viscosity, 
thermal conduction etc.) are bulk properties, whereas in- 
teraction with the environment is a surface phenomenon. 
To this assertion let it be opposed that viscous flow im- 
plies the presence of walls, coupled to the system, where 



the relevant collective momentum is input or withdrawn. 
Furthermore, referring to the celebrated Joule experi- 
ment, where compressed gas is made to expand spon- 
taneously in an evacuated vessel, it is easy to demon- 
strate that the initial acoustic shock requires interven- 
tion of the walls in order to reach final relaxation. The 
softer the walls, the faster does the acoustic perturbation 
vanish. In the unphysical extreme hypothesis where no 
coupling at all to the environment would exist, not even 
that represented by action of the ubiquitous electromag- 
netic radiation field (black body), the system's dynamics 
would be governed by its only conservative Hamiltonian. 
The motion would then retain at all times the memory 
of its initial conditions. 

Considering that the dynamics of irreversible processes 
implies necessarily at some stage action of the surround- 
ings, a general procedure has been developed that quan- 
tifies the latter's most relevant properties. This is con- 
veniently done by referring to thermodynamics, to be 
extended to conditions out of equilibrium. Thermody- 
namics opens indeed the door to the definition of in- 
tensive variables conjugate to every extensive property. 
For given properties, differences of the intensities in and 
out the system measure how much the system and its 
environment are removed from their mutual equilibrium 
conditions. The procedure has been elaborated and dis- 
cussed extensively elsewhere |]llj-|l4|| . It has been applied 
successfully to a manifold of simple relaxing systems and 
to systems supporting steady transport of extensive prop- 
erties (energy, momentum etc.). It has been shown to 
remain valid in conditions far removed from equilibrium, 
even beyond predicted bifurcations. The present work 
expands some of the early results to mixtures of gases 
and compares the theoretical expectations to published 
experimental data. 

In the traditional approach by Boltzmann and his fol- 
lowers, unbound free flow implies the usage of a La- 
grangian description to formalize the motion [[l5|| . That 
is why, using symbol / to represent the distribution func- 
tion in phase space, and omitting extraneous forces for 
simplicity, Boltzmann's dynamic equation is written & 



df 
dt 



df_ 

dx k 



C(f), 



(1) 



where C represents the collision integral. The resulting 
distribution functions are time dependent, even in sta- 
tionary conditions leading to physically less trans- 
parent conclusions. By contrast, by referring relaxing 
systems to their fixed boundaries and considering local 
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thermodynamic properties, the thermodynamic approach 
fits remarkably in an Eulerian frame |l5| , |l7|] 

d i = [m + j. ( 2 ) 

In the latter equation, [/, H] is the Poisson bracket de- 
scribing the implicit motion while J is a source/sink term 
expressing explicit action of the environment. Integrat- 
ing the equation yields very simple expressions for the 
relevant distribution functions. In stationary conditions 
the distribution functions are time independent. This a 
great advantage of the procedure. 

The present contribution is structured as follows. Gen- 
eral principles supporting thermodynamics of systems 
out of equilibrium are outlined in section 2. The for- 
malism leading to prediction of transport properties is 
developed in section 3, where it is applied to single com- 
ponent dilute gases. Extension towards mixtures of the 
atomic gases and comparison with published experimen- 
tal data is the subject of the last section. 

II. THERMODYNAMICS OUT OF 
EQUILIBRIUM 

It is impossible to specify exactly the state of a complex 
macroscopic system (macrostate). We must content our- 
selves with descriptions that are considerably less than 
complete. In fact, our exact information about the prop- 
erties of many-particle systems is restricted to a small 
number of mechanical observables or constraints, directly 
related to the system's Hamiltonian. Among others, let 
us quote: the number of particles of any sort (j): Nj, 
the total energy E and the accessible physical volume V. 
That are the traditional micro-canonical constraints or 
extensive variables. In non- equilibrium conditions, ad- 
ditional constraints prevail, like the total linear momen- 
tum P, the total angular momentum and also possible 
momenta of all the properties cited above. 

A. The Entropy 

Any function determined completely by the set of con- 
straints describing its particular macrostate is a function 
of state. The theoretical definition of the function of state 
entropy goes back to Boltzmann. 

S = k B \n[W(A)}. (3) 

For the inventor, W(-A) meant "wahrscheinlichkeit" 
which is probability. Digging for the realities hidden 
behind this word may lead to some controversies but, 
using the same initial letter, most authors wisely pre- 
fer the English weight of the given observational state 
or macrostate. The latter is interpreted as the mea- 
sure of the domain accessible to the motion in phase 



space, given the set of constraints (represented here by 
the collective variable A) describing the system's particu- 
lar macrostate. An equivalent definition for VF(A) is the 
number of quantum states (= microstate) all compatible 
with the given set of constraints. Let it be stressed that 
Boltzmann's definition is applicable to non-equilibrium 
conditions simply by including the additional constraints 
in the definition. 

Let the list of the extensive constraints defining a given 
macroscopic system in a particular macrostate be written 
{A r }. The entropy is a function of this collection of 
variables. By differentiating with respect to this set we 
get by definition the set of conjugate intensive variables 
or intensities {£,•}■ 

dS 

dS = X! ojr dXr = - kB J2 & dx - ( 4 ) 

r r r 

Equation (4) is Gibbs' celebrated differential equa- 
tion, generalized to possible non-equilibrium conditions. 
It renders the usual temperature (dS/dE) -1 and the 
collection of chemical potentials —T(dS/dNj). In non- 
equilibrium conditions it generalizes the definition by at- 
taching a given intensity to each of the additional non- 
equilibrium constraints. 

Boltzmann's definition of the entropy is valid whatever 
the number of particles in the system of interest. A sig- 
nificant advantage of referring to its differentials, namely 
the intensities, is that their values are independent of 
this number and also of the discrete nature of physical 
systems. 

B. Generalized Massieu function 

If two systems are allowed to exchange some extensive 
properties it is easy to show that the state of mutual 
equilibrium, that is the condition where exchange van- 
ishes on the average, occurs when the conjugate intensi- 
ties equalize |T§| ] . The total entropy becomes then insen- 
sitive to possible infinitesimal fluctuations in the relevant 
exchange [l^ ]. 

Let the two systems to be considered be a huge body 
representing the surroundings (reservoir) on one hand 
and a small object called the system on the other. Their 
respective dimensions are such that thermodynamic flows 
do not alter significantly the reservoir's intensive vari- 
ables. The latter are therefore constants, defining the 
external experimental parameters and fixing the con- 
straints imposed to the smaller system. Intensities are 
indeed better measured and controlled. Therefore, in- 
stead of referring to the entropy, an explicit function of 
the extensive properties ({A,.}), thermodynamics makes 
widely use of thermodynamic potentials and Massieu- 
Planck functions, obtained from the extensive properties 
by Legendre transformations ]ig[] . 

The constant volume Legendre transform of the en- 
tropy is the generalized Massieu function M(V, £ r ). It is 
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obtained by performing the following transformation on 
the entropy: 



S 



M(v,z r ) = —+y2t r x r , 



(5) 



whereby the volume is not included in the collection of 
indexed constraints. 

Differentiating M. with respect to any intensive vari- 
able yields readily the conjugate extensive properties. 



X r = 



dM 



(6) 



C. Dilute Gases 

With ideal gases, the expression for the generalized 
Massieu function takes a very simple form. Individual 
particles assumed to be independent, the global motion 
may be represented by a swarm of points in a reduced 
6iV-dimensional single-particle phase space (T). 

Let /(r) be a particle distribution function. Any ex- 
tensive property X r may then be related to a generating 
function <j) r (T) so that 



X r 



<Mr)/(r)dr. 



(7) 



With that formalism, using the Lagrangian multipliers 
procedure to specify the maximum entropy conditions 
compatible with the set of constraints, it is easy to show 
fl2f that function f(T) becomes 



/(r) = exp [X><Mr)]. 



(8) 



With ideal gases, if equation (5) is implemented with 
the two latter expressions, M. takes the very simple form 



M(tr,V) = jTexpE^(T)]dr. 



(9) 



The numerical value of this function is the (average) num- 
ber of particles contained in the system. Through the 
integration limits in configuration space it has the sys- 
tem's physical dimensions (volume) as one of its inde- 
pendent variables. By restricting the integration to the 
only momentum coordinates, a local generalized Massieu 
function is obtained, the value of which represents the 
average local density in configuration space. 

With real gases, the generalized Massieu function is 
modified due to the interaction potential between the 
particles. The simplified formulation is however still use- 
ful as an approximation in low density conditions, when 
the duration of the inter-particle collisions is negligible 
compared to the time separating collisions. With hard 
spheres this is certainly the case. 



III. TRANSPORT COEFFICIENTS 

One of the main objectives of the theoretical approach 
to non-equilibrium dynamics is the prediction of trans- 
port coefficients from first principles. Comparison be- 
tween the predicted and experimental results is often con- 
sidered as a test for the validity of the relevant attempt. 

Since Boltzmann first proposed his kinetic equation 
there has been a considerable literature concerning the 
calculation of the transport coefficients |2(],|2lj . Most fre- 
quently cited are the traditional Chapman and Enskog 
derivations [ fl6| and the Green-Kubo formalism. 

It has been stressed above that, for all but perhaps 
a few mechanical properties, exchange occurs more or 
less readily with the surroundings, tending to equalize 
the conjugate intensities to the reservoir values. This 
justifies thermodynamic expressions based on intensities. 
When intensities conjugate to exchangeable properties 
are different from the reservoir values, we have transient 
conditions from where the system tends to relax. If the 
system of interest is connected to a surroundings that is 
not at equilibrium, it reaches and remains in a stationary 
state out of equilibrium. This is the condition we shall 
focus on now. 

If the system is interacting with two reservoirs at dif- 
ferent temperatures separated by some distance (here: 
2D), the conditions of the surroundings define and dic- 
tate to the system the genuine non-equilibrium intensity 
"temperature gradient" . Similarly in the Couette flow 
problem, the externally imposed gradient is caused by a 
couple of walls moving in opposite directions. This gen- 
erates in the system a non-equilibrium intensity "gradi- 
ent of shear momentum" . Asymmetric exchange with 
the two reservoirs produces flows. In this section, the 
relevant transport coefficients will be examined using a 
thermodynamic description. For simplicity, the discus- 
sion will however be limited to hard sphere atomic gases. 

In extremely low density systems, where the mean free 
path is comparable or longer than the system's physical 
dimensions (Knudsen gas), properties picked up by any 
particle from one wall are transported in a single jump to 
the opposite wall. Transport is very efficient indeed. In 
the thermodynamic limit (non-Knudsen regime), head- 
on collisions of like particles do not slow down the trans- 
port properties. By contrast, parallactic or off-axis inter- 
particle collisions do. Their effect is one of reducing the 
range of free transport, while information about the con- 
ditions prevailing in the external reservoir and available 
at the boundaries is transferred to the relevant region of 
the bulk. As a result, local values of the thermodynamic 
properties are justified. 

For the same reason, the flow rates depend on the av- 
erage periodicity r of the perturbing collisions. 
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A. Single component gases 

Let us consider an arbitrary property X r with gener- 
ating function </> r (T). We assume that its flow is directed 
along the z-axis. Let us consider a plane positioned at 
coordinate z* . The basic equation for the flow J r of the 
relevant property through this plane is 



0, 



d 3 p 



(z* —pz t /m) 



br(T)exp[J2ZiMr)}dz. (10) 



The symbol d 3 p is a short form for dp x dp y dp z . Integra- 
tion limits in momentum space is — oo — > oo. Plank's con- 
stant in the numerator represents a norm in phase space, 
related to the classic- quantum correspondence fig] . In 
the subsections to follow, this equation will be applied to 
different types of flow. 

It may be useful to stress the difference between equa- 
tion (10) and that proposed elsewhere, in own relevant 
different contexts for flows of extensive properties fuilpil : 



Jr = 1 1 1 m MT) ex PE^( r )]§- 



(11) 



Contrasting with the latter, equation (10) confirms that 
transport occurs during the free motion time separating 
relaxing collisions. The lesser the collision frequency, the 
more effective is the transport. Collisions do not activate 
transport. They do increase resistance to transport. 



1. Viscosity 

We consider a fluid bound by a pair of walls moving in 
opposite directions (Couette flow). Excepting negligible 
higher order corrections (proportional to r 2 (/3mZ3 2 )~ 1 ), 
the system's conditions are defined completely (at the 
lowest order in t) by the set of constraints listed in table 
H The intensities under direct control of the surroundings 
(exchangeable) are the particles number, the kinetic en- 
ergy and the intensity conjugate to the gradient of shear 
momentum. It may indeed be verified readily that the 
velocity of the walls (y- direction) equals ±<7 y //3. Two 
variables remain to be determined, namely 02 and 72. 
(Index 2 refers to quadratic moments). They require two 
independent equations. 

In stationary conditions there is no local accumulation 
of the transverse component of momentum (p z ) (no pres- 
sure gradient, no acoustic perturbation). The relevant 
flow is therefore independent of z*. Likewise, the total 
flow of energy through the system is zero. By imple- 
menting equation (10) with the two relevant generating 
functions, the conditions dJ Pz /dz = and Ju = yield 
together 



5 72 
2 



ma 



(12) 



The flow of shear momentum may now be determined 
by implementing equation (10) with the generating func- 
tion p y , where 62 and 72 have been replaced by their 
values. This yields 



2/3D (3 



(13) 



where n = M /V represents the particle density. 

The phenomenological reaction at the plates' level 
compensating for transfer of momentum from wall to wall 
is friction. Shear viscosity is the ratio of the sum of the 
forces applied to the two plates to the velocity gradient 
(<7 y (/3Z?) _1 ). Following equation (13), its value is 



(14) 



2. Thermal conduction 

Now we consider a system in thermal contact with a 
couple of heat reservoirs at different temperatures sep- 
arated by a distance 2D. The system's stationary non- 
equilibrium conditions are completely described (at the 
lowest order in r) by the set of constraints listed in ta- 
ble |n|. By inspecting the generating function conjugate 
to the temperature gradient it is clear that fc^VT = 
— 7i (/9 2 .D) -1 . (Index 1 refers to linear moments). 

Knowing that the system is bound by a pair of im- 
pervious walls it may look strange that collective motion 
of the particles perpendicularly to the walls needs to be 
anticipated in constructing the expression for the non- 
equilibrium distribution function / (equation 8). 

In a system controlled by two heat reservoirs at dif- 
ferent temperatures, kinetic energy is not uniformly dis- 
tributed among the particles. Those moving towards the 
cold wall have been equilibrated with the system up- 
stream in a hotter region at the instant of their previous 
collision and vice-versa. In moving from the hot wall to 
the cold one, particles travel on the average faster than 
in their return cycle. If the particles are to change their 
average kinetic energy in a correlated fashion on impact 
with either walls, while the container (the pair of walls) is 
to remain immobile, collective momentum is transferred 
by the container into the system. 

The intensities under direct control of the surround- 
ings (conjugate to exchangeable properties or otherwise 
given constant properties) are the intensity conjugate to 
the particle number, the temperature (or better (3) and 
the temperature gradient. Two intensities remain to be 
determined: 0\ and a z . This requires two independent 
equations. One is the condition for stationarity. The 
other equation describes mechanical equilibrium of the 
system between its walls. 
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According to the local description, at any point in the 
bulk of the system, the average particles density n re- 
sponds to the equation 

/>oo roc />oo 

n(z*) = — / dp x I dp y / dp z 

" J —oo J — oo J —oo 

exp{^6.<M(z = z*),p]}. (15) 

r 

At any position z*, we consider the partial density 
n + (z*) of the only particles with positive velocity along 
the z-direction. Stationarity implies that this partial den- 
sity is compensated exactly by the sum of the densities 
of the particles issued from regions from where they will 
be reaching this position undisturbed in one collision pe- 
riod, their velocities being oppositely oriented (p < 0). 
Hence 



1 



dp x I dp y I dp z 

— oo J —oo J —co 

exp{5>«M(z = z*-^),p]}. (16) 



To the lowest order in t, Relation n_ (z*) = n + (z*) yields 



h - 2—) — = 2a z 



(17) 



The second equation expresses position independence 
of flow of momentum across the system. In other words, 
there are no pressure gradients. Equation (10) is used 
with p z as the flow defining generating function. The 
condition dJ p Jdz = yields (to the lowest order in r) 



5 7i 
2/3' 



(18) 



Flow of energy (heat) through the system is given by this 
same general equation (10) where the flow defining gen- 
erating function is now ^(p 2 /2m). For particles associ- 
ated with internal rotational motion (Eucken correction 
jl6],^2]]), the relevant contribution to the energy should 
be added to the latter generating function. With atomic 
gases the result reads 



5 n 
2/?2 



J E = „lw\ a z - 



1(0 - 111)—} 

2 { 1 2 (3> mD 1 



(19) 



By implementing the latter with the relevant values of 9\ 
and er z , flow of energy becomes 



T - 1£ ll T 

E ~ 8 /3 2 D n /3m' 



(20) 



Heat conductivity (A) is defined as the ratio between 
the sum of the rates of heat exchange at either walls 
(2Je) and the temperature gradient. Hence 



A = — k B n — . 

4 [5m 



(21) 



In equations (14) and (21) the transport coefficients are 
expressed in terms of the effective collision periodicity r. 
For direct comparison with experimental results, an ad- 
ditional expression is required that relates the collision 
periodicity to the mechanical properties of the colliding 
species (mass and cross-section) at the given tempera- 
ture. Here there remains an uncertainty concerning the 
model to be adopted for relaxing collisions, yielding the 
effective cross-section as a function of the temperature. 
For evaluating the results presented above it is therefore 
advisable to eliminate the variable r in the discussion. 
That is where the Prandtl number comes in. It is related 
to the ratio between viscosity and heat conductivity and 
given by 



Pr 



yep 

m A ' 



(22) 



with c p as the constant pressure heat capacity. By imple- 
menting this definition with the results obtained above, 
the experimental values are obtained identically, thereby 
corroborating the general model. fllq,f2l|. 



B. Mixtures of atomic gases 

We are now investigating transport processes in binary 
mixtures of gases. Let the components be indexed A and 
B, where A points to the component with the highest 
mass. Each may be considered as a separate system, 
with its own thermodynamic functions, interacting si- 
multaneously with the other component and with the 
environment. The generalized Massieu function being an 
extensive properties, we have for the composite system 



M = M A + Mi 



(23) 



With dilute gases or gases interacting as hard spheres, 
the individual generalized Massieu functions are defined 
as in equation (9). For each component separately the 
generating functions to be used are the same as for single- 
component gases (see tables | and [n]), excepting for the 
additional indexing of the mass of the relevant particles in 
the generating function for kinetic energy. In stationary 
or quasi-stationary conditions (see below), for exchange- 
able properties where equilibrium between the subsys- 
tems prevails, the intensities are the same. In the exam- 
ples treated below, that will be the case for the tempera- 
ture and its moments and for the intensities conjugate to 
the collective motion. Intensities conjugate to the popu- 
lations and their own moments however will be indexed 
according to which component they refer to. 

It has been stressed above that the collision periodicity 
is an essential ingredient in the dynamics of systems out 
of equilibrium. In multi-component systems, there is an 
average collision periodicity for each of the constituents 
(ta>Tb)< It measures for each component how long the 
relevant atoms move freely before being halted by the 
matrix formed by the other particles, making them feel 
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the thermodynamic conditions dictated by the environ- 
ment. 

In multi-component systems there are homogeneous 
and heterogeneous relaxing collisions. Their frequencies 
add up. The efficiency for exchange of momentum from a 
colliding atom to the local thermodynamic bath depends 
on the masses of the collision partners. When a heavy 
particle hits a light constituent of the thermodynamic 
bath, its path is less disturbed and less momentum is 
transferred than in the opposite case. 

We assume a particle with mass mi and linear mo- 
mentum P hitting a stationary matrix particle with mass 
m2. If the exit path of the matrix particle forms an an- 
gle ij) with the incident one, momentum transferred to 
the matrix equals 2|P| cos(ip)ni2 / (mi + 7712). Hence, the 
relative transfer efficiency of heterogeneous collisions is 
2m2/(mi+m2). For the total effective collision frequency 
of atoms of one sort with respect to the matrix (the recip- 
rocal of t), the latter coefficient is the appropriate scaling 
factor relating the efficiency of heterogeneous collisions to 
the homogeneous ones. 

We take the atoms to be hard spheres. The collision 
cross-sections are respectively cIaa, d BB , cIab- Using the 
scaling parameter defined above, omitting the common 
factor jgy/ f3/n and adding for either constituents the 
individual effective collision frequencies, the two total ef- 
fective collision periodicities read 




n B d\ B 



m A + m B 
2mAm B 



, (24) 



mA + m B n B d 
2mAm B 



2 \ -1 

BB 



(25) 



Considering that transport coefficients of mixtures are 
systematically compared or normalized to either one of 
the pure gas values, the common factor cancels. 

For simplifying the formalism it is advisable to replace 
riA by xn and n B by (1 — x)n. This will be done system- 
atically below. 

With hard spheres we have d,A B = (<1aa + d BB )/2. It 
appears that experimental accuracy of the published data 
on the viscosity of mixtures of atomic gases is sufficient 
to allow the heterogeneous hard sphere diameter to be 
corrected by a factor e close to 1. 



1. Viscosity 

We consider a binary mixture of atomic gases bound 
by a pair of walls moving in opposite directions (Couette 
flow). The mole fraction of substance A is written x 
(N A + N B = N). 

The intensities under direct control of the surround- 
ings (exchangeable properties or otherwise independent 
properties) are the intensities conjugate to the particle 
numbers of either substances [a a, &b), the temperature 
(or better (3) and the linear moment of shear velocity 



Oy (see table |) . Three intensities need still to be deter- 
mined, namely the quadratic moment of the temperature 
(or better 72) and the quadratic moments of the particle 
distributions for A and B (02,A> ®2,b)- 

The three additional relations required for completing 
the thermodynamic description of the system are of the 
same vein as those used for Couette flow in single com- 
ponent gases (see above). For symmetry reasons, flow of 
shear momentum is independent of the particular values 
of the additional intensities. The principles involved in 
their determination will therefore be postponed until the 
section concerning thermal conductivity and diffusion. 

Flow of momentum is supported by either components. 
For each, the contribution is given according to equation 
(10), where the generating function to be implemented 
as 4> r is p y . Integration yields 



d i) 



2(3 2 D 



f- [XT A + (1 - X)T B ] 



The viscosity of the mixture is therefore 



Tjn 



[XT A + (1 - x)t B ], 



(26) 



(27) 



where n/f3 is the total pressure (V). 

In figure [j] the result of equation (27), is plotted for 
a mixture of Xe in He. The experimental results at 291 
K published by E. Thornton and coworkers 23 are in- 
dicated on the same graph (experimental uncertainties 
~ ±1%). The correction factor e for heterogeneous col- 
lisions may be estimated by fitting the curve to the ex- 
perimental results. The curve obtained without the cor- 
rection factor (e = 1) is displayed as a dotted curve. 

The same fit has been performed on the ten different 
mixtures of a tomic gases at the same temperature of 291 
K. Table III lists the values of e giving the best result for 
each mixture. 



2. Diffusion and thermal conduction 

We consider now a binary mixture of atomic gases in 
thermal contact with a couple of heat reservoirs at differ- 
ent temperatures separated by a distance 2D. The mole 
fraction of the heaviest substance (A) is written x. 

For each of the two components, the stationary non- 
equilibrium conditions are completely described by the 
set of constraints listed in table 0. The intensities must 
be indexed accordingly. 

The intensities conjugate to the particle numbers of 
either substances (a a, ct B ), the temperature (or better 
0) and its gradient (or better 71 ) are under direct control 
of the surroundings. Thermal interaction between the 
subsystems removes the necessity of indexing the latter 
two intensities. 

Three intensities remain to be determined namely the 
two gradients of the particle distributions (0i,a, 8i,b) 
and the intensity conjugate to the collective momentum 



G 



from wall to wall (<J Z ). Hence, three additional conditions 
or equations are needed in order to describe the system 
completely. 

Two of the additional conditions are identical to those 
discussed for single component systems. It are mechani- 
cal equilibrium and stationarity of the total particle dis- 
tribution. 

Mechanical equilibrium of the system between its walls 
implies vanishing total pressure gradient. It does not re- 
quire per se vanishing pressure gradient for cither sub- 
stances separately. A possible pressure gradient of A is 
neutralized by an opposite gradient for B. The condition 
is formalized by stating that the sum of the contributions 
of either substances to flow of momentum between the 
boundaries is position independent 



dJ, 



dJ, 



Pz,B 



dz 



dz 



= o. 



(28) 



Applying equation (10) with (f> = p z for each of the two 
substances leads to 

x{0i,A ~ 0) + (1 - x)(9 ltB - 0) = 0. (29) 

Let it be stressed that the pressure gradient for A, is 



l.A ~ 77 



5 71, 
2 0- 



(30) 



with ua = xn, and mutatis mutandis for B. 

The condition for stationarity is defined along the same 
lines as above (equations 15-17), where the densities 
n + (z*) and n_(z*) are now understood as the sum of 
the different components. As a result, the relation for in- 
ternal collective motion (er z ) becomes (see equation 17) 
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The last condition to be considered concerns mutual 
diffusion or motion of the subsystems with respect to each 
other. By implementing equation (10) with the generat- 
ing function <p r — 1 the particle flow of either subsystems 
is obtained, according to whether the parameters in the 
exponential function are indexed A or B. The results are 
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The first contribution in either equations (that pro- 
portional to a z ) represents collective drag generated in 
the fluid by correlated effect of the two walls. This acts 
on the two subsystems alike. Therefore it does not drive 



diffusion of one subsystem with respect to the other one. 
By contrast, diffusion is described by the second part of 
the flow equations. As it may be verified, this is driven 
by the relevant partial pressure gradient. 

Diffusion coupled to flow of heat and vice-versa are 
known as the Dufour and the Soret effects Jlq] . Let us 
consider substance A as the solute and B as the solvent. 
If transport of heat and matter are expressed using the 
convenient parameters for the relevant generalized forces 
(the conjugate intensities, here 71 and 0i a), Onsager's 
phenomenological equations are retrieved p4| , p5| . 

The diffusive stationary state is reached by differen- 
tial displacement of the subsystems with respect to each 
other. Then we have for either subsystems vanishing par- 
tial pressure gradients. This represents therefore the re- 
maining constraint for complete thermodynamic descrip- 
tion of the stationary non- equilibrium system. The ques- 
tion is however how fast the diffusive stationary state may 
been reached in practical cases when thermal conductiv- 
ity of multi- component mixtures is measured. 

The experimental procedure for measuring thermal 
conductivity consists in preparing an appropriate bi- 
nary mixture in equilibrium conditions in a conventional 
thermostat, the mixture being then introduced between 
two walls held or brought at different temperatures. 
When mechanical equilibration is established (relaxation 
of acoustic perturbations), the total pressure is flat (equa- 
tion 29). Nevertheless, on transferring the mixture in the 
region with the temperature gradient, individual pressure 
gradients on the subsystems may have been created, forc- 
ing the particles to segregate. If the mixture consists of 
particles with different mobility it is expected that estab- 
lishment of final stationary conditions is slow. The slower 
moving particles tend indeed to remain distributed ho- 
mogeneously, as they were before the establishment of 
the temperature gradient, while the partial pressure gra- 
dient of the faster moving subsystem compensates for the 
former's resulting partial pressure unbalance. 

When comparing predicted values of thermal conduc- 
tion to experimental data, there is uncertainty as to 
how much the system has been allowed to relax the 
slow coupled particle segregation in the relevant mea- 
surement. Let us assume this would not have occurred at 
all (pseudo-stationary state). The two subsystems may 
then be considered as acting independently for all the 
properties concerning the particle distributions. They re- 
main however tightly coupled for all the properties that 
are promptly interchanged. In particular, they share the 
same value of and 71. The intensity a z conjugate to the 
collective momentum generated by the temperature gra- 
dient is also common to the two subsystems. Concerning 
the latter, its relation to the other intensities and to the 
collision periodicities is given by equation (17). Instead 
of equation (31) we have now two relations, namely 
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By combining the equation for overall mechanical sta- 
bility (equation 29) with the two latter ones, an expres- 
sion for the gradients of the individual partial pressures 
may be derived. Writing 
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this relation reads 
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In practical cases, when heat conductivity is measured, 
the system may be somewhere between the two extreme 
conditions. The uncertainty concerning how close diffu- 
sion has reached stationarity in the experimental condi- 
tions where the measurements have been performed ex- 
plains why thermal conduction data of mixtures are dif- 
ficult to reproduce. Let us express the uncertainty by a 
coefficient c to multiply the right-hand side of equation 
(37). When discussing a homogeneous set of data with 
varying compositions x, we assume for simplicity that the 
same coefficient is valid. 

Transport of heat is supported by either components 
of the mixture. 
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external control of macroscopic systems. It may there- 
fore be claimed that the procedure presented above is 
characterized by a greater physical transparancy when 
compared with the traditional treatments |lq 1 . The num- 
ber of adjustable parameters it implies is low and their 
physical meaning is straightforward. 

The theory confirms that inter-particle collisions are 
not responsible for the transport processes. By contrast, 
transport occurs during the free motion time of the parti- 
cles in the periods separating collisions. By colliding, the 
particles exchange their mechanical properties with the 
matrix or bath of the remaining particles representing 
the system, the local thermodynamic properties of which 
are defined by the external constraints. Permanent inter- 
action of this matrix with the environment causes corre- 
lations between individual motions to disappear. The 
effect of collisions is to increase the resistance opposed 
by the system to flows. 

The theory relies on the hard sphere mutual interaction 
model. This implies further the definition of an effective 
temperature dependent collision diameter for the differ- 
ent particles involved. The theory is clearly not suitable 
for predicting this quantity. It requires trajectory calcu- 
lations based on the particular inter-particle interaction 
pot entials. The question is well documented elsewhere 
P6J. In the treatment presented above, this problem has 
been bypassed by focussing on results obtained at con- 
stant temperature. 



For each, the contribution is given according to equa- 
tion (19), where the relevant intensities are as determined 
above. Hence, 



1 



7 7i \ T A 



2 (3 ttiaD 



(1 — x)n 

0* 



7 7i x t b 



2 (3 ttlbD 



(39) 



In comparing the result with experimental data, coef- 
ficient c may be taken as an adjustable parameter. This 
exercise has been performed on the data at 291 K pub- 
lished by E. Thornton and coworkers p6|. The coeffi- 



cients yielding the best fit are listed in table IV. Figure 
^ is an illustration of the results. The accuracy is better 
than the announced experimental precision (4%). 



IV. CONCLUSION 

The theory for transport of extensive properties in di- 
lute mixtures of atomic gases has been successfully devel- 
oped in the context of the thermodynamic approach to 
non-equilibrium processes |l2j ]. The equations are based 
on accurate definitions of intensities conjugate to the set 
of extensive properties defining the system's particular 
macrostate. Intensities are the natural parameters for 
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TABLE II. List of the main constraints for thermal con- 
ductivity (distance between the walls: 2D) 



X r 


0r(r) 




Particles number 


1 


a 


Gradient of particle distribution 


(z/D) 


01 


Kinetic energy 


XX 2m! ) 


-P 


Gradient of energy distribution 




-7i 


Collective transverse momentum 


Pz 





TABLE III. Correction e to the heterogeneous collision di- 
ameter of pairs of gases, obtained by fitting the viscocity of 
the relevant mixtures to experimental data (temperature: 291 
K) 





Ne 


Ar 


Kr 


Xe 


He 


1.03 


1.01 


1.00 


0.97 


Ne 




0.98 


0.98 


0.98 


Ar 






1.03 


1.00 


Kr 








0.99 



FIG. 1. Predicted viscosity of a mixture of Xe in He, with 
e = 0.98 (smooth curve) and published experimental data 
(temperature: 291 K). The dotted curve is for e = 1. 



FIG. 2. Predicted thermal conductivity of a mixture of Xe 
in He taking c = 0.5 (smooth curve) and published experi- 
mental data (temperature: 291 K). The dotted curves are for 
c — and c = 1 



TABLE IV. Correction c to the partial pressure of the pairs 
of gases, caused by unrelaxed thermal diffusion, obtained by 
fitting the predicted termal conductivity to experimental data 
(temperature: 291 K) 





Ne 


Ar 


Kr 


Xe 


He 





0.7 


0.55 


0.5 


Ne 




0.5 


1 


1 


Ar 






1 


1 


Kr 








1 



TABLE I. List of the main constraints for Couette flow 
(distance between the walls: 2D) 

X r 4>r(r) £r 

Particles number 1 a 

2nd moment of particle distribution [(z/D) 2 — 1] 62 

Kinetic energy XX 2k) ~P 

2nd moment of energy distribution [(-§) 2 — 1] XX 2^7) — "f' 2 

Gradient of shear momentum (z/D)p y a y 
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